Nanoparticle size distribution estimation by full-pattern powder diffraction analysis. 
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The increasing scientific and technological interest in nanoparticles has raised the need for fast, 
efficient and precise characterization techniques. Powder diffraction is a very efficient experimental 
method, as it is straightforward and non-destructive. However, its use for extracting information 
regarding very small particles brings some common crystallographic approximations to and beyond 
their limits of validity. Powder pattern diffraction calculation methods are critically discussed, with 
special focus on spherical particles with log-normal distribution, with the target of determining size 
distribution parameters. A 20-nm Ce02 sample is analyzed as example. 

PACS numbers: 61.10.Nz, 81.07.Bc, 61.46. +w 



I. INTRODUCTION 

We are assisting at a booming expansion of nanoparti- 
cle research and technology. Synthesis method especially 
make fast progresses^. Analysis methods, however, are 
not up to speed. A fundamental simple task as deter- 
mining and controlling the size distribution of nanoparti- 
cles (NPs hereafter) is currently a complex experimental 
work, involving electron microscopy and combined tech- 
niques. In this work we want to highlight the possibilities 
offered in this issue by a much less complex technique as 
powder diffraction. 

Powder diffraction is a widespread technique with a 
great potential to meet the increasing demands of mi- 
crostructural material characterization. The methods 
of powder diffraction data analysis have reached matu- 
rity for micrometer-sized polycrystallinc materials. How- 
ever, when the particle size falls much below 100 nm, 
specifically tuned methods of analysis are needed to ex- 
tract meaningful information from powder diffraction 
patterns. In fact, nanoparticles (NPs hereafter) present 
unique analytical challenges. In the most complex cases, 
non-crystallographic structures^i^SiLSiSiiSiii may oc- 
cur. Surface-related deformation fieldsiSiiii^ are an- 
other challenge. In these extreme cases, the classical 
crystallographic formalism becomes quite useless. The 
Debye scattering function^ (that is, the direct evalua- 
tion of the NP structure factor from the interatomic dis- 
tances) is the only choice in those cases. We are currently 
developingifl*i& methods to increase the efficiency of such 
calculations and make them a practical tool. 

Even for crystalline NPs, however, the small size plays 
a decisive role. Bragg peaks may be so much broad- 
ened that they cannot be simply separated and many 
approximations, commonly accepted for micrometer size 
domains, fail. As we will show, also models specifically 
corrected for NPslii&i 9 . may fail for ultra-small NPs (say 
below 5 nm diameter, as it will be better specified). 
Again for these ultra-small sizes the Debye scattering 
function is the only choice for obtaining precise results, 



while the smaller number of atoms makes it extremely 
practical. 

The plan of the paper is the following. In Sec.[n]we dis- 
cuss the shape-based method for calculating NP powder 
patterns in relation to the surface structure and to its lim- 
its of validity at small sizes. Application to full-pattern 
fit on a test-case (20-nm CeC^) is shown in Sec. IIIII 42 
Summary and conclusions are given in Sec. II VI 



II. POWDER PATTERNS AND SIZE 
INFORMATION 

Schcrrer's formula™ is the most known method for ex- 
tracting size information from powder patterns (namely, 
from the Bragg peaks' width). This is a simple method, 
but accurate only to the order of magnitude. However, 
since Scherrer's work, line profile analysis has made enor- 
mous progresa 21 ! 22 ! 23 ! 24 ! 25 ! 26 : 27 ! 28 ! 29 ! 30 . 

Theoretical progress on understanding the physical ori- 
gin of peak broadening has been focused on the disloca- 
tion analysis, size broadening being considered as a side 
effect to be corrected for in order to determine the de- 
fect structure. Nevertheless, today it is possible to deter- 
mine the parameters of a (log-normal) size distribution of 
crystallites, together with information on type and con- 
centration of dislocations. These methods are, however, 
complex and sophisticated, requiring a fairly high signal- 
to-noise ratio, low and flat background, a precise decon- 
volution of the instrumental broadening and especially 
well-isolated Bragg peaks. 

Full-pattern fitting methods (c/. Sec. Ill All arc more 
direct and robust, especially when the target is the size 
analysis. Firstly, they use all the experimental informa- 
tion, regardless of partial or total peak overlap, increasing 
redundancy and therefore precision and decreasing exper- 
imental requirement. Furthermore, they allow the eval- 
uation of a NP-characteristic feature, namely the varia- 
tion with size of the lattice paramete»i2iil (an effect that 
can be important below 20 nm). Corrections for texture, 
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microabsorption, anisotropic elastic peak shifts and in- 
strumental broadening can also be implemented. 

An efficient and precise method to evaluate NP diffrac- 
tion patterns is needed to perform full-pattern fits. Here- 
after we discuss the shape-based methodiiiSiiS with a 
thorough analysis of its validity limits. 

A. NP shape-based diffraction models 

We shortly recall some methods for the calculation of 
the powder diffraction intensity for a NP with known 
periodic structure and definite size and shape. In the 
following the length of a vector v will be denoted by 
v. Accordingly, q will be the scattering vector of length 
q = 2 sin 0/ A, where 9 is the scattering half-angle and 
A the incident wavelength; h shall denote the scattering 
vector associated with a Bragg peak, its length being 
h. A NP occupies a geometrical region of space G. We 
recali2i*2£ the definition of a shape function S(r), such 
that S(r) = 1 if r lies inside G, S(r) = otherwise. We 
shall hereforth suppose that S(—r) = S(r) so that its 
Fourier transform is real. 

However, defining the shape of a crystal means also 
to describe what happens to the atoms on the surface. 
These are increasingly important at very small sizes. In 
fact, there are different ways of interpreting the action of 
S(r), the most meaningful ones being: 

a) truncating sharply the scattering density (the elec- 
tron density for x-rays) at the surfac e) 31 ' 32 ; 

b) selecting all whole unit cells whose origins arc in G 
and all whole atoms whose centres lie in the selected 
cellsS; 

c) selecting all whole atoms whose centres are in G. 

Useful illustrations are found in Fig. I of Ref. (see 
Figs, fa, Ic and Id, respectively for a, b, c)4£ To evaluate 
the diffracted intensities, in cases b), c), one may utilize 
the Debye function. In this way the chosen model is 
faithfully represented. It is possible, however, to proceed 
in a different way, that is, by the shape-function method. 
Accordingly, we first evaluate the scattering amplitude 
A(q). The explicit expressions^ are, for cases a,b,c: 

A a (q) = Y S(q-h)F(h), (1) 

heA* 

A\q) = Y, S(q-h)F(q), (2) 
A c (q) = J2 S(q-h)F(h,q), (3) 

heA* 

where A* is the reciprocal lattice; S(q) is the Fourier 
transform^ of S(r), or 

S(q) = I 6?rS{r)c 2mq r = f d a rc 2mq r , (4) 



and it satisfies S(q) — S(—q) because S(—r) = S(r); 
F(h) is the unit cell structure factor 

F(h) = YU(h)* 2vih ' ra . (5) 

a=l 

where the sum index a runs on the atoms in the unit cell, 
which have form factors f a and position vectors (relative 
to the cell origin) r a ; F(q) is the same as the former but 
evaluated in q; and F(h, q) is the mixed expression 

N a 

F(h,q) = Y / U<l)e 27rih - r «- (6) 

a=l 

It is evident that form a) is simpler but by construction 
less reasonable - for electron and x-ray diffraction - than 
b) and c). In fact, the sharp truncation of the electron 
density at the surface is unjustified. For neutron nuclear 
elastic scattering the atoms are point scatterers, there- 
fore, construction a) coincides with c). Accordingly, in 
the neutron case, the atomic form factors are constant 
and A a {q) = A c (q). 

Form b) depends on an appropriate choice of the unit 
cell. Clearly, it preserves the stoichiometric composition 
and symmetry. 

Form c) needs a careful implementation (regarding the 
definition of G) to preserve stoichiometry, that is im- 
portant for ionic compounds; however, it is clearly more 
flexible. Remark also that, in the case of monoatomic 
lattices, instead - as for simple-cubic, face-centered or 
body-centered cubic metals - construction b) and c) will 
be coincident and A h (q) = A c (q). 

B. NP scattering intensities 

Squaring Eqs. HIL'l.'il) we obtain the intensities. Sup- 
posing S centrosymmctric and S real, we have 

I a {q) = J2 S 2 {q-h)\F{h)\\ (7) 

heA* 

I b (q) - \F(q)\ 2 Y, S 2 (q-h) 7 (8) 

F(q) = Y S 2 (q-h)\F(h,q)\ 2 . (9) 

heA* 

Here, we have neglected cross-summations of the form 
TZ(q)= Y S(q-h)S(q-k)M^ h Ml k (10) 

h,fe£A* 

where ovcrbar stands for complex conjugate and, for 
x=a,b,c, respectively, it is M* h = F(h), M% h = F(q) 
or Mg h = F(h,q). Neglecting IZ(q) is, first of all, a 
question of convenience, because its evaluation - cither 
analytical or numerical - is a nightmare. 
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There are obvious reasons for neglecting IZ(q) for large 
particles. Consider a spherical particle with cubic struc- 
ture with lattice parameter a and radius R 3> a. S(q) 
is large only for q < 1/R, and decreases as {2nqR)~ 2 for 
q^> 1/R. As for any Bragg peak h it is 1/R <C 1/a < h, 
K(q) ~ 0((R/a)- 2 ) can be neglected. 

Fo r smal ler particles the situation is different. In 
Refs. Il7l33l it is proposed that IZ(q) is negligible due to 
a certain statistical 'smearing' of the NP surface region 
on a thickness of the order of the lattice parameter a. 
However, this hypothesis cannot be accepted by default. 

Firstly the order at the surface strongly depends on 
the considered crystal phase and on the actual sample. 
Consider that for a NP of diameter D = Na, the frac- 
tion of atoms included in a layer of thickness a is w 6/jV 
(about 50% at D = 10a, still 12% at D = 50a). The 
structure of this large fraction should be carefully consid- 
ered on a case-by-case basis. Relaxations in the core due 
to a disordered layer of thickness a should also be con- 
sidered. Secondly supposing a default smearing of the 
NP boundaries flattens the different construction princi- 
ples of forms a, b, c. In fact, the differences among them 
regard the finest details of the NP surface structure. 

We shall hereafter assess the effect of neglecting 7Z(q) 
on the calculation of a powder diffraction pattern. In 
App. B we carry out some relevant calculations. Evidently 
this will depend on the choice of form a, b, or c. Examples 
are reported in the following section. 

For form I h (q) it turns out that, even when 7Z(q) is not 
negligible, it yields a contribution that is approximately 
proportional to the retained term I h (q) of the scattered 
intensity. This means that the effect of neglecting 7Z(q) 
may be just a small error on the global scale factor for 
samples composed of particles of equal size. However, as 
this effect is size-dependent, it may hamper the evalua- 
tion of size distribution when this is not very narrow. A 
size-related correction factor for the scale factor may - 
and should - be evaluated (see App. J in this case. This 
of course is an undesired complication. 

In cases a) and c) the neglected term IZ(q) depends on 
the crystal structure (see App.^l. It is not a constant scale 
factor for all Bragg peaks, and it may have a significant 
gradient in the Bragg peak positions. At very small sizes 
the latter may induce a systematic error also in the lattice 
constant determination. However, in the x-ray case, for 
form a) IZ(q) is larger - and has a larger gradient in the 
Bragg peak neighbourhood - than the corresponding term 
for form c). 



C. NP powder patterns 

To obtain a powder diffraction pattern, we must 
integrate I x (q) (x=a,b,c, see Eqs. illl-?i:>l) ) at con- 
stant q. We write q in polar coordinates as q = 
(g sin t/j cos </>, g sin -0 sin </i>, gcos-0) = (g, w), where a; is the 
orientation defined by the pair (ip,4>)- We have to inte- 
grate over the set of all orientations £1 = {0 < tp < ir, 0< 



4> < 2tt} (with duj = sin ipdipdcf)) , as 



sini/>d</> / d0/*(g) = g 2 / dw/ x (g,w). (11) 
Jo Jn 

In detail, considering the expressions for the different 
cases, we have 



TO = q 2 



£ \F(hf f du>S 2 (q-h); (12) 
= <1 2 £ / du\F(q)\ 2 S 2 (q-h); (13) 



i c M) = q 2 



£ l^,?)| 2 / duS 2 (q-h). (14) 

h£A* ^ n 



The integration in case b) is much more difficult and it 
cannot generally be expressed in closed form even for 
simple shapes. Therefore, as a careful implementation of 
form c) is at least as good a description as form b), we 
shall disregard b) in the following. Suppose now that G 
is a sphere of radius R and volume V = 4ttR 3 /3, we have 



S(q) = S(q) = W 



sm (y) ~ ycosjy) 

y 3 



(15) 



y=2nqR 



(16) 



and, as \q — h\ = (q 2 + h 2 — 2qhcosip) 1 / 2 

s(q -h) = w \ Mv)-v<*»M 

y 

with y = 2ir{q 2 + h 2 - 2qhcosip) 1/2 R. 
Substituting in Eqs. (|12I14|I yields^ 

W - nsr £ |F(, *> |2 (jL t M -- (17) 

heA* 

w - inM)| 2 ^H(18) 

where A± = y~ 2 (l - sin(2y)/y + sin 2 (y)/y 2 ) 
for y = 2irR(q±h). 

Now we consider the crystal's Laue group Q so that we 
can extend the summation on the asymmetric part A* JQ 
of the reciprocal lattice: 



i c P (q) = 



3qVR V11 \F(hM 2 {A -~ A+) 

heA*/g 

^— ^fi h \F(h,q)\ 



(19) 



(20) 



hGA*/g 



where [ih is the multiplicity of h subject to Q . Evaluation 
of Ip(q) is only slightly more complex than I~{q), and the 
gain in accuracy justifies the effort. 

We have computed test patterns to compare forms a) 
and c), considering NPs of diameter ps 10a, being this the 
lower size limit of validity of the shape-based approach. 
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We have considered Au spherical NPs of diame- 
ter 5 nm (a=0.40786 nm, A=0. 154056 nm, 26 = 
20° . . . 150°, Lorentz correction and Dcbye- Waller factor 
cxp(-5g 2 /2), with B = 0.005 nm 2 ). The powder pat- 
tern was calculated exactly by the Debye sumiS^ and 
by Eqs. (|19I20J) . The profiles showed in Fig. GJi are cal- 
culated on an absolute scale. They match quite well, but 
a maximum error «2 - 3% is present in both cases a,c. 
The profile wR agreement index between /^ )cbyc and 7° is 
3.1%, between /^byc and J£ is wR=AA%. The difference 
profiles (Fig.^>) show that I® chyc — Ip has a similar shape 
to J° ebye , while J° ebye - I* is quite different. Accord- 
ingly, refining a scale factor between /Dcbye anc j jc j owers 
wR to 2.0% (with featureless difference, Fig. [Tp), while 
a scale factor between /Dcbyc anc j ja yi e ^ s uiR=3.5%, 
with still a characteristic difference profile. Furthermore, 
the peak positions result very little shifted (< 0.002°) 
between J° ebye and 1°, while they are shifted up to 0.04° 
between i"° ebye and ja (Fig.[TJi). 



Then, we have considered ZnSe spherical NPs of di 
ameter 4.8 nm (a=0.5633 nm, A=0. 154056 nm, 26 
20° . . . 135°, Lorentz correction and Debye- Waller factor 
with B = 0.005 nm 2 ). Once more, the powder pattern 
was calculated exactly by the Debye sumi2ii£ and by 
Eqs. i|19l20|) . The profiles - calculated on an absoultc 
scale (Fig. |2Jt) - match quite well with a maximum error 
~ 1 — 2% for both cases a,c. The profile agreement in- 
dex wR between 1^° and J£ is 1.8%, between I^obyc 
and is wR=3.l%. The difference profiles (Fig. 0d) 
show again that I® ehye — 7° has a similar shape to I® chyc , 
while /Dcbyc _ ja j g q U j^ e different. Accordingly, we have 
refined again a scale factor (and this time also a differ- 
ent Debye- Waller factor B) between I^ chyc and 7°. wR 
decreases to 1.6% with featureless difference (Fig. |2p). 
On the opposite, when refining scale factor and Debye- 
Wallcr factor between /^ cbyc and I* the agreement in- 
dex does not go below wR=3.1%. Also the difference 
profile is little changed (Fig. |2p) . Again, the peak posi- 
tions result very little shifted (< 0.001°) between I° ebye 
and n, while peak shifts up to 0.05° between /Dobyo an ^ 
Ip are visible (Fig. |2Ji) ■ Form c) again turns out to be 
less affected than a) by neglecting the cross-term 1Z. A 
small variation of the Debye- Waller factor (from 0.005 to 
0.0047 nm 2 ) is due to the fact that the 7?.-neglection error 
changes slightly the intensity ratios. This is however less 
troublesome than the peak shifts observed for form a). 



It results that at NP diameters D sw 10a the errors in 
the shape-based diffraction pattern calculations, what- 
ever form we choose, start to be evident. This approach 
should not be used below this threshold. Also, form a) - 
which is the standard choice for large particles - shows a 
much larger error and should be avoided in favor of c). 



D. The log- normal size distribution 



There are several experimental and theoretical 
reasonsM to believe that NP powders have a log-normal 
distribution of NP size. The log-normal distribution of 
NP radii is usually written in terms of its mode R m and 
width Wr, as 



L(R) = 



1 



w R V2n 



exp 



[l0g(g) - log(/i T1 

2w% 



(21) 



The most direct information on a distribution is provided 
by the distribution-averaged NP radius R aV e and the rel- 
evant standard deviation <tr. For a log-normal, the latter 
parameters are related to the former by 



R a 

and 



J R m exp( W 2 s /2); a R = R 2 m exp(2w 2 ); (22) 



Rn 



y/1 + al/Rl, 



l0g(l + 4/^e). 



_ We shall use a form depending directly on R a 
Setting two adimensional parameters p = R/R a 



(23) 
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1 + a R /R 2 ave , we have 
L(R) = 



iV 27rl °S (c) 



exp 



log 2 (pV~c) 



21og(c) 



(24) 



Volume- and area-averaged NP diameters can be derived 
by 



3 o 
- R r 3 - 



D 



±R c 2 

^ 1 h ave^ 



(25) 



III. NANOCRYSTALLINE CERIA 



A. Experimental 



X-ray powder diffraction 
nanocrystallinc 20-nm Ce02 
able for a round-robi: 



.35 



patterns of a 
sample, avail- 
were downloaded 



(http : //www. boulder .nist .gov/div853/balzar/, 
http://www.du.edu/~balzar/s-s_rr.htm). The NP 

size is well inside the limits of validity of the shape-based 
method. Among the available datasets, the selected raw 
data were collected at the NSLS X3B1 beamline of the 
Brookhaven National Laboratory in flat-plate geometry, 
with a double-crystal Si(lll) monochromator on the 
incident beam (A = 0.6998 A, 26 = 12°(0.01°)60°) and 
a Ge(lll) analyzer crystal on the diffracted beam. 



B. Data preprocessing 

Three data preprocessing stages have been accom- 
plished. First, the instrumental function has been de- 
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convoluted by an original advanced technique, includ- 
ing denoising and background subtraction, described in 
Ref. HE 

Secondly, the pattern has been fitted by generic asym- 
metric Voigt profiles so as to obtain information about 
peak positions and intensities. By comparing the inten- 
sities as evaluated from the fit with the theoretical ones a 
small correction for texture and/or microabsorption has 
been evaluated. The intensity corrections so obtained 
have then been stored and used in the subsequent stages. 

Finally, the peak positions were found to be slightly 
anisotropically shifted. This has been attributed to a 
small residual stress, due e.g. to dislocations. To con- 
firm this point, we have evaluated the average lattice 
spacing variations [Ad/d]h = — jfo cot(8f l )A(29h) for all 
single reflections h. Then we have compared those val- 
ues with a simple model of elastic anisotropy2£. They 
resulted in good agreement. In Fig. E3 we show the fit of 
[Ad/djh's with Eq. (28) of Ref. HI The magnitudes of 
the residual stress tensor components, at least for those 
which can be determined in this way, resulted to be in 
the range 1-10 MPa. The values of |A(20 h )| are below 
0.005°, and [Ad/d]h range in l-7xl0~ 4 , which are quite 
small values. As the strain broadening is of the same 
order of magnitude of the peak shifts^, we can confirm 
that strain broadening is rather small in the Ce02 sample 
and can be neglected, as in Ref. |3f|. Also the residual- 
stress peak shifts so obtained have been saved as fixed 
corrections for the subsequent stages. 

C. Full- pattern refinement 

The total intensity diffracted by the powder NP sample 
is described by the sum 

k ma x 

jcai {q) = Ibkg{q) + g L {R k )I k {q), (26) 

k=0 

where R k = (k + l/2)AR, k = l...k max ; I k (q) is 
I*{q) of Eq. J2DJ} evaluated at R = R k ; and I bkg (q) 
is a polynomial modelling the background. The step 
AR = a(27r/3)~ 1 / 3 is chosen so as to have an integer 
number of atoms in each fc-th x-ray sphere of radius 
R k , while keeping the point density constant and pre- 
serving stoichiomctry. It is evidently possible to use a 
size-dependent lattice parameter a k in the calculation of 
Ik (<?) ■ For this sample this has been deemed unnecessary. 
Indeed, for diameters of 20 nm, the lattice parameter of 
CeC>2 has been founds to be already equal to the bulk 
value. A least-square full-pattern refinement means min- 
imizing the quantity 

N obs 

x 2 = J2(l cal (q i )-lf s ) 2 w i . (27) 

i=X 



Here 1° is the i-th point of the experimental pattern 
corresponding to the scattering vector qi, N t s the num- 
ber of experimental points and the weights Wi are the 
estimated inverse variance of the observations. The re- 
fined parameters are: the average NPs radius (R a ve) and 
the radius dispersion or, the isotropic Debye- Waller fac- 
tors B for O and Ce atoms, the cubic unit cell parameter 
a and seven background coefficients. For the minimiza- 
tion, we have used (for this work) a modified simplex 
algorithm^, which is robust but time-consuming; how- 
ever, computing times were reasonable. A derivative- 
based algorithm (Newton, in progress) should give a 
handsome acceleration. 

The final results are given in Tab. U together with the 
corresponding values of Ref. The Debye- Waller fac- 
tors result to B Ce = 0.0065 nm 2 and B = 0.0084 nm 2 . 
The calculated profile is plotted in Fig. 21 with the ex- 
perimental pattern and the profile difference. The excel- 
lent fit quality and the final GoF value (1.21) indicate 
the achievement of a reliable result. Indeed, the esti- 
mated parameters are in good agreement with Ref. l35l 
The slight discrepancy (sa 0.2 nm), larger than standard 
deviations, might be explained by the improved deconvo- 
lution method here applied and by the use of the whole 
pattern instead of a limited number of peaks as in Ref.l35l 



TABLE I: Comparison of size distribution results. Standard 
deviations are in brackets. Units are nm. 





This work Ref. jfr 


Rave 


9.58 (0.02) 9.33 (0.07) 




4.138 (0.003) 3.92 


a R.I Rave 


0.1866 (0.0008) 0.177 0.003 


D v 


24.01 (0.02) 22.8 (0.4) 


D A 


17.98 (0.02) 17.2 (0.2) 



IV. CONCLUSIONS 

The method of shape-convolution to calculate the 
diffraction pattern of NP powders has been thoroughly 
discussed with respect to its limits of validity. Concerns 
in applying this method below its optimal size range have 
been demonstrated theoretically and by simulated pat- 
terns. Finally, the effectiveness of full-pattern powder 
data analysis based on the shape-convolution method 
was proved to obtain precise size distribution information 
on NP powder samples with a log-normal distribution of 
spherical crystallites. 
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Figures 



FIG. 1: (Color online) A model spherical cluster AU3925 of 
5.0 nm diameter with fee structure (a=0. 40786 nm) has been 
constructed according to principle c) of Sec. Ill Al In this 
CeLS6 ; £LS the monoatomic fee Wigner-Seitz unit cell contains 
one atom, principle c) coincides with b). 

a: the powder diffraction pattern: red, exact intensity 
Ip Bhye (q) calculated by the Debye function^; blue dotted, 
Ip(q) calculated by approach a), Eq. 119H : green dashed, 
Ip(q) calculation by approach c), Eq. 12011 . All intensities have 
been calculated on an absolute scale and then scaled by the 
same factor. 

b: lower line, red, difference Ip ye (g) — Ip(q); middle line, 
green, difference Ip^ ebyc (q) — Ip(q); upper line, blue, the exact 
powder pattern x 1/100 (Debye method) for comparison, 
c: lower line, red, difference I p e (q) — s Ip(q) after refin- 
ing a scale factor s = 1.024; middle line, green, difference 
Ip ehye (q) - sIp(q)aheT refining a scale factor s = 1.027; upper 
line, blue, the exact powder pattern x 1/100 (Debye method) 
for comparison. Note that the c)-type pattern difference is 
flattened while the a)-type retains sharp contributions, 
d: detail around the (111) peak of the 7° ebye and Ip patterns 
(same coding as in part a) after scaling, showing a significant 
peak shift for the Ip pattern. 
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FIG. 2: (Color online) A model spherical cluster (ZnSe)i289 of 
4.8 nm diameter with fee structure (a=0.5633 nm) has been 
constructed according to principle c) of Sec. lH Al In this case, 
as the fee Wigner-Seitz unit cell contains two atoms, construc- 
tion c) differs from b). 

a: the powder diffraction pattern: red, exact intensity 
7° ebye (g) calculated by the Debye function^; blue dotted, 
Ip(q) calculated by approach a), Eq. 11911 : green dashed, 
7p(g)calculation by approach c), Eq. 1201 . Again all inten- 
sities have been calculated on an absoulte scale, 
b: lower line, red, difference 7^ )cbyc (g) — Ip(q); middle line, 
green, difference 7^ >ebye (g) — Ip{q)\ upper line, blue, the exact 
powder pattern x 1/100 (Debye method) for comparison, 
c: lower line, red, difference 7° ebye (q) — s Ip(q) after refining 
a scale factor s — 0.999 and an overall-isotropic Debye- Waller 
factor B = 0.468 (7° obyc has been evaluated with B = 0.5); 
middle line, green, difference 7° ebyc (g) - sl^(q) after refining 
a scale factor s = 0.996 and B = 0.467; upper line, blue, the 
exact powder pattern x 1/100 (Debye method) for compari- 
son. Note again that the c)-type pattern difference is flattened 
while the a)-type retains sharp contributions, 
d: detail around the (531) peak of the 7° obyc and Ip 1 patterns 
(same coding as in part a) after scaling, showing a significant 
peak shift for the Ip pattern. 
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FIG. 3: (Color online) The measured (red error bars) and 
calculated (blue diamonds) lattice spacing variations for all 
well-isolated Bragg peaks of Ce02 plotted against the rele- 
vant peaks diffraction angle. Error bars have been evaluated 
assuming a constant error of 0.0006° on the anisotropic angu- 
lar peak shift. Calculated values refer to the model of Ref. 1371 
where residual stress components have been refined. 
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FIG. 4: (Color online) Nanosized Ce02 powder pattern final 
fit. Blue diamonds - the observed deconvoluted intensity; red 
continuous line - the calculated intensity; black continuous 
line, below - difference profile (same scale). 
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APPENDIX: ERROR EVALUATION 

Assume to deal with particles of centrosymmetric 
shape S(r) = S(—r) and equivalent spherical radius R 
(i.e., the radius of the sphere of equal volume.). The 
shape Fourier transform S(q) is then a real even func- 
tion: 

§(q) = S(-q). (A.l) 
Recall also that the gradient of an even function is odd: 



G(q) 



V q S(q) 



-G(-q). 



(A.2) 



Our aim is to evaluate - for the different forms a), b), 
c) as introduced in Sec. Ill Al and carried out in Sec. Ill Bl 
Sec. Ill CI - the neglected residual intensity contribution 
TZ(q) of Eq. i|10|) with respect to the respective retained 
term (cf. Eqs. Q7I8I9[0 in the immediate vicinity of a 
Bragg peak. 

Let ho the nearest Bragg peak to q. First note that, if 
\q — h \ ^> 1/R, TZ(q) is of order (2nqR)~ 4: , so we neglect 
it altogether. If q is very close to ho, set q = ho + Aq 
(so Aq < 1 /R) . We can drop in the sum over h all terms 
with h ^ h because they are 0(2i:qR)^ i and reorder 
the second sum, obtaining 



TZ a {ho + Aq) » S(Aq) S{Aq + k)F{h ) 



x F(h a - fe) 



K h (ho + Aq) 



S(Aq)\F(h + Aq)f 
S(Aq + k) 



(A.3) 



(A.4) 



fcGA* 
fc^O 



U c (ho + Aq) « S(Aq) J2 S(Aq + k) 

feeA* 

X F(h Q ,q)F(h -k,q) 
where q = \h + Aq\ 



(A.5) 



At the same time, for q = ho + Aq with Aq < 1/R, 
the intensities I x (q) of Eqs. (|7I8I9J) can be approximated 
by the ho-th term of the RHS sum, neglecting terms of 
0(2TrqRy 4 . Furthermore, in general, S(Aq) = 5(0) + 
0(Aq 2 ). Therefore, the ratios lZ K (h + Aq)/r(h Q + Aq) 
are given by 



n a (h ,Aq) = 



K a (h Q + Aq) 



S 2 (Aq)\F(h )\ 2 



s(A9 ' fc)4 pfer (A - 6) 



fc£A* 
fc#0 



n h (h ,Aq) 



TZ c (h ,Aq) = 



K b (ho + Aq) 
= S2(Aq)\F(ho + Aq)\2 
* £ §{Aq,k) 

fc^O 

TZ c (ho + Aq) 
S 2 (Aq)\F(h ,q)\* 

Y ,?(Aq,fc)5^M 
^ 1 ' F(h ,q) 



(A.7) 



(A- 



fc^O 



where 5(Aq, fc) 



5(Aq + fc) 

5(0) 



Note that, because of Eqs. (|A.1IA.2|I . we have 

S(Aq,k) = s(-Aq,-fe); (A.9) 
q(Aq,fc) = V Aq s(Aq,k) = G(Aq + k)/S(0) 

= -g(-Aq,-k). (A.10) 

We can immediately veryfy that in case b) it results 

W Aq TZ b (h Q ,Aq) = 9(Aq,k) (A.ll) 

feeA* 
fc^o 

In the sum above the term with index k is always accom- 
panied by a term with index — k. Setting also Aq = 0, 

and using Eq. jmnji, we have 

V A<J £ b (/i o ,0) = (fl(0,fc)+fl(0,-fc)) = 0. 

fceA*/2 

(A.12) 

where A*/2 denotes an arbitrarily chosen half-space of 
the reciprocal lattice without the origin. Now, expanding 
7Z h (ho, Aq) in Taylor series at Aq = 0, we have 



K h (ho,Aq)^K h (ho,Q) + 0(Aq 2 



(A.13) 



Note also in Eq. (|A.7f) that TZ b (h o ,0) does not depend 
on the considered Bragg reflection ho- Therefore, we can 
write 



TZ h (q) (X I h (q), 



(A.14) 



and the proportionality constant can be evaluated by 
Eq. fP|) with Aq = 0. We can conclude that the ef- 
fect of neglecting lZ h will be just a relative error on the 
global profile scale factor. This factor is size-dependent, 
however, therefore for size distribution analysis at small 
sizes it may be necessary to introduce a correction as 
from Eq. (|A.7|I . 

Cases a), c), are more complex. We are interested 
to powder diffraction, where I(q) is to be integrated at 
constant q, therefore we shall consider 



K X (h ,Aq) 



1 r 
2 



K x (ho,Aq) +K X 



-ho, —Aq) 

(A.15) 
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Expanding TZ (ho,Aq) in Taylor series at Aq = 0, we 
have 

It (h , Aq) « K*(h Q , 0) + W Aq Tt(h , 0) • Aq + 0(Ag 2 ). 

(A.16) 

We shall now develop TZ (h ,Aq) and V & q 1Z (h n ,0) in 
cases a, c. 



Then, evaluating the gradient in Aq = 0, using Eq. (|A.2|) . 
we have finally 

V Aq ft>o,0) = 9(0, k)x 

fceA*/2 

x Re[r(/i ,fe)-r(h ,-fe)](A.24) 



1. Case c) 

First, recall that the atomic form factors f a (q) are con- 
stants for neutron scatering and monotonically decreas- 
ing smooth functions in the x-ray case. In the latter 
case, furthermore, the form factors of different elements 
have remarkably similar profiles. For a structure with N a 
atoms in the unit cell, it is then possible^ to approximate 



1 AT. 

f a (q) « c a (f(q)) = c a — 2J fp(q) 



(A.17) 



13=1 



with c Q appropriate constants. Therefore the structure 
factor ratios appearing in Eq. i|A.8|) can be simplified as 



F(h - k,g) _ J2 a =i c^e 
F(h ,q) 



-27ri{ho—k)-r 



2irihn -r 



o-r l3 



independent of q = \ho + Aq\. Note that 



= T(h ,k), 

(A.18) 



r(-h , -k) = r(h ,k); r(-ho,k) = r(h , -k). 

— ( A - 19 ) 
Now we can write explicitly TZ using Eqs. (|A.8IA.15p 

and 

lT(h ,Aq) = 1 [s(&q,k)T(h ,k)+ 

fc^O 

+s(-Aq,k)T(h ,-k)}. (A.20) 

Splitting the sum, reordering k — > — fc in one part, using 
Eq. (|A.9J| and recombining, we have 



TZ C (h ,Aq)= s(Aq,fe)Re[r(/i ,fe)] ■ (A.21) 

fceA* 

Again as in Eq . Eg) , we can pair terms with fc and 
— fc. Using Eq. (|A.19|) . we obtain 

7T(h ,Aq) = ^ Re[s(Aq,fc)T(h, ,fc) + 

feeA»/2 

+«(Ag,-fc)T(/io,-fc)]. (A.22) 

Define now the arbitrary half-lattice A*/2 as that defined 
by a plane _L ho passing through the origin and contain- 
ing ho. The origin is excluded. We have 

7f(/i ,0)= s(0,fe)Re[r(/i ,fe)+r(h 0) -fe)]. 

fceA*/2 

(A.23) 



The gradient V AqTi (ho,0) is a vector. We have to take 
its angular average to determine the effect on the powder 
pattern. This is done by simply taking the scalar product 
with ho = ho/ho'- 

^A q Tt(h Ql O)-ho = (9(0,k)-h )x 
fceA*/2 

x Re[r(/i ,fe)-r(h ,-fejjA.25) 



For spherical shape, it will be g(Q, fc)||fc; therefore 
terms with fc _L ho will be zero and those with fc||fao 
will be most important. Both g(Q, fc) and s(0, fc) 
arc damped oscillatory functions with amplitude ~ 
(2nkR)~ 2 . As 1/a < fc, the magnitudes of both TZ (h ,0) 
and V &qlZ (ho,0) are of order (a/R) 2 . Unfortunately, 
Eq. I|A.25|) cannot be estimated more in detail, because 
of the dependence from the 'reduced' structure factors 
r(ho,k). However, we can assess that its importance 
would be smaller than the corresponding term for case 
a) for x-ray scattering. 



2. Case a) 

In case a), we can trace the same steps as in case c) 
but instead of the 'reduced' structure factors r(ho,k) we 
have to consider the ratios 



«ho,k)= F tr k) = F{h °:y5 ho) (a.26) 



F(h Q ) 



\F(h )\< 



and in the analog sums of Eq. ljA~23)l and Eq. (|A~25)) for 

TZ (ft-0,0) and VaqT?. (^o,0) there will appear terms as 



Re [£(/*>, fc)±C(/»o,-fc)]. 



(A.27) 



The most important terms for the powder pattern are 
again those with fc||/io- The structure factors F(ho ± fc) 
(see Eq. (J3J) depend on form factors f a (\ho±k\), and for 
k\\ho these will be strongly different. This in turn will 
amplify the differences £(ho,k). Therefore it is likely 
that for case a) the effect of the neglected term TZ a will 
be significantly larger than for case c). The examples 
reported in Sec. Ill Cl show iust that. 



